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Is Cr02 Fully Spin-Polarized? - Analysis of Andreev Spectra and Excess Current 
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We report an extensive theoretical analysis of point-contact Andreev reflection data available in 
literature on ferromagnetic CrCh- We find that the spectra can be well understood within a model 
of fully spin-polarized bands in Cr02 together with spin active scattering at the contact. This is in 
contrast to analyses of the data within extended Blonder-Tinkham-Klapwijk models, which lead to 
a spin polarization varying between 50 % and 100 % depending on the transparency of the interface. 
We propose to utilize both the temperature dependence of the spectra and the excess current at 
voltages above the gap to resolve the spin-polarization in CrC>2 in a new generation of experiments. 

PACS numbers: 74.55.+v,74.45.+c,72.25.Mk 
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Half-metallic fcrromagncts arc materials with one of 
the two spin bands metallic and the other insulating. 
Their characterization has attracted great attention, 
since a fully spin polarized ferromagnetic material can be 
very useful for fabricating spin batteries and ideal mag- 
netic tunnel junctions used in spintronics applications |l( . 
There are only a few materials that are suspected half- 
metals [2|, one of them being CrC>2. 

Following early experiments by Soulen et al. [3[ and 
by Upadhyay et al. [J], point-contact Andreev reflec- 
tion (PCAR) has been extensively used to probe the 
spin-polarization of strong fcrromagncts. In this method 
a nano-sized point contact is formed by pressing a su- 
perconducting tip into the ferromagnetic material. The 
conductance- voltage characteristics is recorded and com- 
pared with theory in order to extract the polarization of 
the ferromagnet. The key ingredient in these experiments 
is the suppression with increased spin polarization of An- 
dreev reflection processes. Andreev reflection is the scat- 
tering event at which an electron quasiparticle incident 
from a non-superconducting metal is retro-reflected as a 
hole quasiparticle in the opposite-spin band. Charge con- 
servation is upheld by injection of a Cooper pair into the 
superconductor. With increased spin-polarization such 
retro-reflection is suppressed. 

In recent experiments [a, Iff , supercurrent was observed 
to flow in long Josephson junctions of &O2. If Cr02 is 
half-metallic, these observations indicate that a conver- 
sion of supercurrent carried by spin singlet Cooper pairs 
in the superconductors to supercurrent carried by triplet 
(equal spin) Cooper pairs in the halfmetal is taking place. 
We have [7| proposed a conversion mechanism based on 
spin-active interface scattering. The question arises, if 
this model stands the test of other experiments. Here, 
we show that spin-active interface scattering can explain 
also the PCAR experiments [3|, Isl— ll 51] on &O2 within a 
model that assumes fully spin-polarized bands. 



The conventional analysis of the PCAR experiments 
(see Fig.[T]for the setup) relies upon an extended Blonder- 
Tinkham-Klapwijk (BTK) formula [16j, consisting of two 
terms. The first term includes Andreev scattering and 
is the usual BTK formula [17| for a point contact with 
an unpolarized material. The second term is the con- 
ductance of a point contact with a completely spin po- 
larized material for which Andreev reflection is absent. 
The two terms are weighted according to the formula 
G = (1 — P)Gn + PGh, where P is the transport spin 
polarization. Taking the limits V —¥ and T — > 0, 
the conductance would reach G(0) = 2(1 — P)G n for an 
ideal contact without backscattering (unit transparency) 
since then Gn — > 2G„ while Gh = in the whole sub- 
gap range. Here, G n is the conductance in the non- 
superconducting state. The BTK model includes an 
interface barrier modelled by a delta-function potential 
quantified by a dimcnsionless parameter Z, where Z = 
corresponds to unit transparency and Z 3> 1 to low trans- 
parency. The fit of the conductance-vs-voltagc data re- 
sults in values for the barrier strength Z and the polariza- 
tion P. It has often been necessary to also use the super- 
conducting gap A as fit parameter. Since the gap is a fit 
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FIG. 1. Point-contact between a superconductor (SC) and 
ferromagnetic Cr02 , including a spin-active scattering region. 
The arrows indicate a contact effective magnetic moment and 
the bulk magnetization. 
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TABLE I. Results of non-linear curve-fits of the extended BTK model (middle set of columns) and the spin-active interface 
model (right set of columns) to point-contact Andreev reflection data on CrC>2 with superconducting tips of Nb (A(T = 0) = 1.5 
meV, T c = 9.2 K) and Pb (A(T = 0) = 1.35 meV, T c = 7.2 K). We have obtained improved fits to the extended BTK model 
[lq ] by including a series resistance r s — R s /R n normalized to the point contact resistance R n as a fourth fit parameter [14[ [. 
The resulting spin polarization P, barrier strength parameter Z, and zero temperature gap parameter A are therefor different 
than found in the original papers. In the spin-active interface model the bulk quantities are fixed (P=100 % and A retains its 
bulk value), while two interface parameters (barrier strength Z and spin mixing angle i?) and the series resistance r s have been 
used as fit parameters. The report in Ref. [15| only contains the fit parameters, but no spectra. 



parameter, the gap feature in the spectrum is shifted in 
an unpredictable way. A shifted peak position has been 
attributed to either a suppressed superconducting order 
parameter near the contact (shift to lower voltage) or a 
spread resistance (shift to higher voltage). The spread 
resistance R s is the resistance of the material between 
the contact (the tip) and the voltage probe which should 
have been eliminated if a true four-point measurement 
could have been set up [14|, see Fig. Q] 

It is very unsatisfactory that the existing analyses of 
the PCAR experiments [aS-El with the extended BTK 
formula have given a wide spectrum of spin polarizations 
in CrC>2 ranging from 50 % to 100 %. In most experi- 
ments the barrier strength is low, with a Z between and 
2. This sample to sample variation leads to a most likely 
spurious dependence of the polarization (a bulk mate- 
rial property) on Z (an interface property). It has been 
argued that the intrinsic polarization can be obtained 
by extrapolating the dependence P(Z) to Z = (for 
example Ref. [ill), although this procedure is question- 
able since the functional dependence P(Z) is unknown 
|l4( . The suppression of spin polarization from the ex- 
pected P = 100 % has in some cases been attributed to 
the unknown interface region. Here, we provide a model 
for this effect. Another unsatisfactory feature of the ex- 
tended BTK model is that the gap A must be used as fit 
parameter. Indeed, well characterized superconducting 
STM tips display the bulk gap [l8|. 

In a different experimental set-up [19|, Zeeman split 
conductance curves of CrC^/Al junctions with fabricated 
tunnel barriers were measured with the Meservey-Tcdrow 
technique [20(. The observed simple linear shift of the 



spectra with applied field strongly points to complete 
spin polarization in &O2, in sharp contrast to the pic- 
ture emerging from the analysis of the PCAR data with 
the extended BTK theory. 

We give here an alternative interpretation of the PCAR 
experiments on Cr02, where we at the outset assume a 
spin polarization of 100 % at the Fermi level and then fit 
the data with two parameters: besides Z we utilize the 
spin-mixing angle $ that describes the difference in scat- 
tering phase factors picked up by electrons of opposite 
spin. This spin-mixing leads to a mixture of singlet and 
triplet Cooper pairs in the superconductor close to the 
contact. If spin-flip scattering is present at the contact, 
triplet correlations are induced also in Cr02 [2J. Spin- 
flips can be induced for example by having misaligned 
moments at the contact, as depicted by the light blue 
contact area in Fig. [T] where the magnetization is mis- 
aligned with respect to the bulk &O2 2l|. When dis- 



cussing PCAR, this physics is closely related to a spin-flip 
Andreev reflection process that leads to enhanced sub- 
gap conductance, here simply tuned by the fit-parameter 
i9. The formula for the conductance relevant for a half- 
metallic ferromagnet with a spin-active interface has been 
derived within the quasiclassical theory of superconduc- 
tivity and is presented as Eqs. (153)-(157) in Ref. [22 1 . 

We have made a fit to all available PCAR data on 
CrC>2 with the spin-active interface model |23[. To make 
the comparison between the two scenarios consistent, we 
have also made new improved fits of the data to the ex- 
tended BTK model of Mazin et al. [16[ , since a variety of 
extended BTK models have been used in the literature. 
The resulting fit parameters are presented in Table 1 and 
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FIG. 2. Point-contact Andreev reflection data on Cr02 
from the literature and non-linear curve fits to the theory of 
superconductor-halfmetallic ferromagnet point contact with 
spin-active interface (red solid lines) and the extended BTK 
model (blue dashed lines). The blue dotted lines are fits to 
the modified BTK model without using A as fit parameter. 
Panels (a)-(d) are the data sets 1, 2, 7, and 10 in Table 1. 



four representative data fits are shown in Fig. [2] 

In Fig. [2Ja) we show the data of Soulen et al. [3| to- 
gether with the excellent fit to the spin-active interface 
model (red solid line). The blue dotted line is the fit to 
the extended BTK model without using the gap as fit pa- 
rameter. Clearly, the fit is not good. On the other hand, 
by letting A vary freely the fit can be improved except 
for about 10 data points in the low-voltage region (blue 
dashed line, A reduced by 33 %). The overall fit, the x 2 
measure, is satisfactory in this case. In FigJ5Jb) we show 
the fit to a set of data from DeSisto et al. [8[ . In this case 
the modified BTK fit can be made excellent by reducing 
A by 10 %. Both models give satisfactory x 2 measures, 
although three data points at low voltage are not per- 
fectly fitted by the spin-active model. In Figs, [^c)-(d) 
we show two other representative fits to both models. 
Again, in (d) the gap had to be reduced by about 30 % 
to improve the fit to the modified BTK model. 

From the obtained fit parameters presented in Tabic HI 
we can draw several conclusions. The main finding of this 
paper is that from a non-linear curve fit (x 2 ) perspective, 
all data sets can be fit well with the spin-active interface 
model with 100 % spin polarization at the Fermi level 
of Cr02- This provides a complementary picture to the 
one emerging from the traditional extended BTK model 
fit. This is also reflected in the same order of magnitudes 
in the calculated average variances a = x 2 /-^ (where N 
is the number of data points in each set) for the two 
models, shown in Fig.EJd). The spin-mixing angle, here 
quantifying spin-flip Andreev reflection at the interface, 
is varying between samples but is suppressed in high- 
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FIG. 3. Correlations between fit parameters, (a) Spin polar- 
ization P and series resistance r s versus the barrier strength 
Z in the modified BTK model, (b) The superconducting gap 
A versus Z in the modified BTK model, (c) Spin mixing an- 
gle $ and r s versus Z in the spin-active interface model, (d) 
Comparison of the fitting function \ 2 divided by the number 
of data points N for the two models. 



Z junctions. This is consistent with the picture of the 
origin of the spin-mixing effect at interfaces to strong 
ferromagnets [7|, where the minority spin electrons gain 
an additional phase compared with majority spins during 
reflection inside the classically forbidden region in the 
half-metallic ferromagnet. As the barrier between the 
materials is enhanced, this effect is reduced. 

The wide spread of the spin polarization P resulting 
from the extended BTK model (50 % to 100 %) is unsat- 
isfactory, since P is a bulk property. We find a correla- 
tion between spin-polarization P and the barrier strength 
Z, see Figj3ja), where (as reported before) P appears 
to approach 100 % as Z — > 0. We have utilized the 
(normalized to the contact resistance) spread resistance 
r s = R s / R n as fit parameter in both models, although it 
plays a more important role for the fit with the extended 
BTK model. The variation of r s with barrier strength 
is shown in Fig. OUa) and (c). Only in one reference has 
the spread resistance been discussed seriously (Woods et 
al. [14j|). However, because we found an uncerta inty 
concerning their data we do not present a fit here [24J ]. 
It is crucial to use A as fit parameter in the extended 
BTK model, while such variation does not improve the 
fits to the spin-active interface model. It is unclear why 
the (bulk) A should vary so much [see Fig. [3Jb)], in con- 
trast to other experiments with superconducting STM 
tips [18(. All point contacts appear to be highly trans- 
parent with a small barrier strength parameter Z < 1. 
This is consistent with a Fermi velocity mismatch, which 
for Pb/Cr02 was estimated in Ref. [l5j to give Z m 0.26. 
But this phcnomcnological parameter should include a 
range of effects causing mismatch between the materials. 
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FIG. 4. Temperature dependence of the zero-voltage conduc- 
tance as predicted by (a) the modified BTK model and (b) 
the spin-active interface model. Excess current at V S> A/e 
versus barrier strength parameter Z as predicted by (c) the 
modified BTK model and (d) the spin-active interface model. 



mentally so far is the excess current, formally defined as 
lexc = liniy^oo [I(V) — I n (V)], where I n (V) is the cur- 
rent in the normal state (= R n V according to Ohm's 
law). In certain limits, the excess current can be com- 
puted analytically although the formulas are rather cum- 
bersome. We present the excess currents predicted by the 
two models in Fig.[4tc)-(d). A measurement of the excess 
current in addition to the PCAR spectrum can be used 
to pin down one of the fit parameters (or as consistency 
check) in future experiments. 

In conclusion, a number of PCAR spectra of Cr02 have 
been presented in the literature where, by comparing the 
data to extended BTK models, a putative spin polariza- 
tion between 50 % and 100 % has been extracted. This 
is in contrast to Zeeman split conductance measurements 
where 100 % polarization was found. We have provided 
an alternative view of the PCAR data, where the spin 
polarization is 100 %, but the scattering at the contact 
is spin-active. 



We are able to fit the experimental data with three 
parameters describing the interface (d and Z) or the ge- 
ometry (r s ), while the extended BTK model relies for the 
same spectra on four fit parameters among which two (P 
and A) pertains to bulk properties, one to the interface 
(Z), and one to the geometry (r s ). The part of the spec- 
tra hardest to fit to either model is the low- voltage region. 
In this region there are typically much less data points 
than in the high voltage region [for example FiglJJb) and 
(c)]. This typically happens in a current bias set-up, 
which is not ideal for PCAR. Our fits could be further 
improved if we would allow for broadening in the form 
of a convolution with a Gaussian (as used in Ref. [Iff), 
which would describe e.g. voltage fluctuations. Since we 
do not know all experimental uncertainties we leave this 
question open for future experiments. 

We would like to point out a few details of importance 
for future PCAR experiments, which may shed more light 
on the properties of Cr02- Thermal smearing is impor- 
tant, since it gives a considerable increase of G(V = 0) as 
compared with the T — >• limit at the temperatures used 
in the experiments. For the spin-active interface model, 
G(0) — > as T — >• independently of the barrier strength. 
This is a unique feature that has not been fully explored 
experimentally. In our model this is a result of vanishing 
spectral current j s = at the Fermi energy e — 0. In 
contrast, in the extended BTK model, G(0) saturates at 
a value given by the polarization and barrier strength, 
see Fig. HJa)-(b). Thus, the temperature dependence of 
G(0) in a well-defined voltage bias set-up can be used as 
a consistency check between experiment and theory. 

Another quantity that has not been explored experi- 
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Our model includes a variation of the spin-mixing angle 
with impact angle 9 P in accordance with a 5-function 
potential of strength Z: $( d v ) — $(0) cos 2 (6 p ). 
The data shown in Ref. [1J| , Fig. 4, is presented as being 
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suming a Pb counter electrode. 



